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In this work we confront dark matter models to constraints that may be derived from radio 
synchrotron radiation from the Galaxy, taking into account the astrophysical uncertainties and 
we compare these to bounds set by accelerator and complementary indirect dark matter searches. 
Specifically we apply our analysis to three popular particle physics models. First, a generic effective 
operator approach, in which case we set bounds on the corresponding mass scale, and then, two 
specific UV completions, the Z' and Higgs portals. We show that for many candidates, the radio 
synchrotron limits are competitive with the other searches, and could even give the strongest con- 
straints (as of today) with some reasonable assumptions regarding the astrophysical uncertainties. 



I. INTRODUCTION 

Dark Matter (DM) is one of the most important issues 
in particle physics and cosmology and understanding its 
nature will likely play an essential role in our compre- 
hension of both fundamental interactions and the struc- 
ture of the Universe. Over the years, a Weakly Inter- 
acting Massive Particle (WIMP) has emerged as one of 
the favourite candidates, in part due to the natural ex- 
planation of its cosmological abundance through thermal 
freeze-out (see e.g. [IH3])- While, as of today, we have 
no other evidence for DM than through its gravitational 
manifestations, the alleged weak interaction of WIMPs 
have open the possibility to actually observe DM experi- 
mentally. Several strategies have been proposed, and are 
actively pursued, to search for WIMPs. Direct detection 
experiments, such as CDMS g] and XENON100 0, are 
dedicated to the search of DM in the vicinity of the so- 
lar system. These are supplemented by multi-purpose 
particle physics experiments at colliders, most notably 
the LHC, where DM is expected to be produced, and to 
manifest itself as missing energy, in collisions. In partic- 
ular, analysis of single-photon or mono-jet events with 
missing energy have recently proved to give very perti- 
nent constraints on WIMP mass and interactions [6rfl0] . 
A radically different, and a complementary approach is 
to search for indirect detection of WIMPs, through the 
remnants of their annihilation (or decays) in astrophys- 
ical environments, like the Galactic Centre (GC) of the 
Milky Way, nearby dwarf spheroidal galaxies (dSphs) or 
in general any dense region of the Universe. The possible 
remnants, or messengers, are high energy neutrinos, anti- 
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matter in cosmic rays (CR), and gamma-rays, or more 
generally an injection of energy of charged particles in 
the early universe. 

In [TT] some of us analysed bounds on effective cou- 
plings between dark matter and SM particles from single- 
photon and mono-jet signals, at LEP and the LHC re- 
spectively. Such studies are usually compared to exclu- 
sion limits set by direct detection experiments (see also 
e.g. [2 El HO]). In general terms, the colliders data 
are comparatively more constraining for low mass dark 
matter candidates, in particular below the threshold of 
direct detection experiments, while the latter are more 
constraining at higher masses, where dark matter pro- 
duction at colliders is impeded. Interestingly, indirect 
searches tend also to be most constraining for low mass 
dark matter candidates. This is essentially because the 
flux of particles produced by dark matter annihilation 
(gamma-rays, etc.) is proportional to the inverse DM 
mass to the square. In particular, interesting constraints 
have been set on the annihilation cross section of DM 
based on the measured synchrotron radiation from the 
inner regions of the Milky Way. While synchrotron radi- 
ation constraints on DM have already been much studied 
in the literature (e.g. in [I2HT6] ). to our knowledge no 
analysis of specific particle physics models implications 
have been made so far. Concretely in the present work 
we confront the constraints from synchrotron radiation 
in the Galaxy at radio frequencies to those set by collid- 
ers data. In particular we study to which extend they 
are complementary. Our analysis is based both on an 
effective operator approach (so we put limits on energy 
scales) and on two specific DM models, the so-called Z' 
and Higgs portals. We also show how colliders and radio 
synchrotron radiation limits compare to bounds set by 
Fermi LAT based on dSphs |17j . and to constraints im- 
posed on DM annihilations from the effect on the CMB 
anisotropies [IB] , 

The constraints from colliders and those from indirect 
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searches do not quite stand on the same footing. In par- 
ticular, although the radio data are potentially strongly 
constraining for rather light DM candidates - as we show 
also in this work - the modelling of DM induced radio 
fluxes suffers from several sources of astrophysical uncer- 
tainties. We explore those by using both a semi-analytic 
approach, which allows to control, for instance, the de- 
pendence of the radio flux on the magnitude of the mag- 
netic fields in our region of interest, and a full numerical 
calculation as implemented in the GALPRDFQcode [191120] 
which allows to explore the full set of CR propagation pa- 
rameters and up-to-date energy losses, thus cross check- 
ing the semi-analytic method and calibrating its param- 
eters. To gauge the impact of CR propagation param- 
eters on synchrotron signals, we sample a range of CR 
propagation parameters sets using both those derived to 
probe uncertainty in the local CR fluxes, and used tradi- 
tionally to bracket this type of uncertainty (MIN, MED, 
MAX set of parameters [21]), and CR propagation sets 
which were recently shown to give a good description of 
the gamma-ray [22l [23] and radio 24] data and which are 
therefore suited to test electron propagation parameters 
in the inner regions of our galaxy. 

We present our results as a function of the systematic 
uncertainties on the modelling of the backgrounds due to 
standard astrophysics, illustrating potential future im- 
provements of the limits based on radio data. 

The paper is organised as follows. In section [II] we re- 
view the formalism of the semi-analytical approach used 
in this work, and also present the setup for the full nu- 
merical study. In Section |III| we discuss at lengths the 



astrophysical framework. In part III A we discuss the as- 



trophysical uncertainties, while in part |IIIB| we discuss 
our cross-check between semi-analytical and numerical 
approaches. In part |III C] we show the dependence of the 
synchrotron flux with different choices of Dark Matter 
profile, and Part |III D| is devoted to study, numerically, 
the impact of different choices of CR diffusion and prop- 
agation parameters to the synchrotron signal. Part |III E] 
is a semi-analytical (cross-checked by numerical) study of 
the dependence of the synchrotron flux on the magnetic 
field normalisation. Then, section [V] is devoted to the 
discussion of synchrotron constraints on particle physics 
models. In part VA we adopt a generic, also known 



as model-independent, approach based on effective op- 
erators. There we confront the synchrotron constraints 
we derive with other indirect searches and colliders con- 
straints. In parts VB and V C we consider two specific 



UV completions, the Higgs and Z' portal, to which we ap- 
ply the synchrotron constraints. While the conclusions 
bear some resemblance, the constraints are quite distinct 
for the two models. Finally we give our conclusion and 
prospects in|VI[ 



II. SYNCHROTRON RADIO EMISSION IN THE 
MILKY WAY 



In the very high frequency radio (VHF) to microwave 
frequency range (i.e. 10 MHz — 300 GHz), various as- 
trophysical processes contribute to the observed diffuse 
emission. The radio sky at frequencies below ~ 20 GHz 
is dominated by synchrotron emission of CR electrons 
(accelerated in e.g. supernovae shocks) on Galactic mag- 
netic fields. However thermal bremsstrahlung (free-free 
emission) of electrons on the galactic ionised gas also con- 
tributes in this range. At higher frequencies the Cosmic 
Microwave Background (CMB) and its anisotropies rep- 
resent the main signal. Emission by small grains of vi- 
brating or spinning dust becomes relevant at even higher 
frequencies, starting approximately at > 60 GHz. 

Synchrotron signal of electron by-products of dark 
matter WIMP self-annihilation is generally expected to 
fall in the 100 MHz —100 GHz range (see for example 
[14] ). However, it has been noted in [T5] [TB] that for 
light dark matter candidates it is beneficial to consider 
lower frequencies. For instance, for a magnetic field of 
fiG strength, as typical for the Galaxy, electrons with 
energies < 10 GeV generate synchrotron emission which 
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peaks at tens of MHz (see Figure 4 in [TB] or Eq. ( 13 ) 
below), and following [16] we focus on a 45 MHz survey 



Additional motivation to consider the 45 MHz survey 
is that it covers a large region of the sky (for a visual 
representation of the coverage of this and other radio sur- 
veys, see [H]). In the past, observations of the Galactic 
Centre compact radio source (SgrA*) with high angular 
resolution surveys have been used to constrain a possi- 
ble DM signal (but at higher frequencies, for example a 
408 MHz survey [Ej] was used in [HQ]!] and a 330 MHz 
observation of the SgrA* [28] in [I!]). Indeed the syn- 
chrotron flux is expected to be strongly enhanced close to 
the Galactic centre, as cusps in a DM density are gener- 
ically found in DM simulations in central regions of DM 
halos, and therefore strong constraints on DM annihi- 
lation cross section may be placed by considering radio 
emission from SgrA*. However, the mass in the inner re- 
gions of the Galaxy is baryon dominated and it is possible 
that baryonic feedback processes might erase the putative 
dark matter cusp [TH3]. On top of that, the very pres- 
ence of a black hole might disrupt the DM profile, both 
by making it steeper, through adiabatic contraction [29] 
or cored through scattering of dark matter particles by 
stars in the dense stellar cusp observed in the vicinity of 
the SgrA* [3U]. Having these uncertainties in mind, we 
take advantage of the fact that 45 MHz [25] is a large 
scale survey and we test its data for a DM contribution 
farther away from the Galactic Centre, where the DM 
density is more robustly determined. 

In order to compute the synchrotron signal, the propa- 
gation and energy losses of a Galactic electron population 
need to be modelled. 
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The propagation of electrons^] in the galactic medium 
is governed by a transport equation, which can be written 
as [3T] 



dn(x, E) 

at 



q(x, E) + V ■ [K xx Vn(x, E) - V c n(x, E)] 



h p2Kpp h^ n{x ' E) 

d_ 

dp 



(pn(x,E)-^(VV c ) 



(1) 



where n(x, E) is the number density of electron per unit 
energy, and q(x, E) is the electron source term. The 
transport through magnetic turbulence can be described 
by the diffusion coefficient K xx . In the following we will 
assume, as customary in literature, that the diffusion co- 
efficient is spatially independent and has an energy de- 
pendence of the form K xx = K E S . V c is the convection 
velocity and re-acceleration is described as diffusion in 
momentum space and is determined by the coefficient 
K pp , and is usually expressed using Alfven speed v a , de- 
fined as 



K K 



3J(4-5 2 )(4-(5) 



(2) 



Cosmic rays propagate in the diffusive halo which is 
usually approximated to have a cylindrical shape with 
radius Rh around 20 kpc and half-thickness Lh which 
could lie in the range of I to 15 kpc. The spatial 
boundary conditions assume free particle escape, i.e. 
n{R h ,z,p) = n(r,±L h ,p) = 0. 

We dedicate the next two sections to describe a semi- 
analytical and, briefly, fully numerical approach used to 
derive synchrotron flux. While semi-analytical approach 
allows in some cases for more physical insight (for in- 
stance, we use it to find the analytical dependence of 
synchrotron flux on the strength of the magnetic field, see 
III E I , we have checked the robustness of that approach 



against a full numerical computation, as implemented in 
the GALPROP code [151120]. 



A. Semi-analytical approach 

For the semi-analytical calculations of the synchrotron 
signal, following [3U [33J we neglect re-acceleration and 
convection terms, which we have checked it is a safe as- 
sumption for our particular studjj^] and rewrite Eq. |l]) 



2 DM annihilation produce the same amount of electrons and 
positrons, so we refer collectively to these particles as 'electrons'. 

3 However, these two parameters are by default included in the 
numerical calculation, see the next section. 



dn(x, E) 
dt 



- V • [K(x,E) Vn{x,E)} 
= q{x,E) , 



dE 



\b{x,E)n{x,E)\ 
(3) 



where b(x, E) encodes the energy loss rate. Cosmic ray 
electrons loose energy mainly through synchrotron radi- 
ation and Inverse Compton scattering (IC), with a rate 
b(x, E) which at the galactic medium is typically of the 
order 10 -16 GeV-s -1 . Additional bremsstrahhmg losses 
of electron energies on the interstellar medium are ne- 
glected in this semi-analytical approach, although they 
may have a very small effect (see for example Fig[I]), 
since for the synchrotron frequency we work with, the 
electron energies of interest are around 1 GeV. 

Assuming a steady state, Eq. ([3| can be re-expressed 



dn(x, E) 
dt 



- K Q L\n{x,E) = q{x,E) 



(4) 



where the variation with energy has been parametrized 
in terms of the parameter t = — f dE(E s /b(x, E)). If at 
the level of propagation one considers that energy losses 
have an average value over all the diffusion reg ion, i.e. 
b{xTE) « b(E) = E 2 /t GeV-s" 1 (see Eq. Q), then 
in Q h(x,E) = E 2 n{ Xl E) and q{x,E) = E 2 ~^q{x,E). 
The solution of this equation can be found in the Green 
function formalism to be 



n(x, E) 



I 



dE s 



HE) j E 

[ d 3 xG(x,E <- x s ,E s )q(x s ,E s ) 
Jdz 



(5) 



where the volume integral is over the diffusion zone (DZ). 
The Green function G(x, E <— x s , E s ) gives the probabil- 
ity for an electron injected at x s with energy E s to reach 
x with energy E < E s , and has a general solution of the 
form 



G(x,E <- x s ,E s ) = -j^G(x,t<- x s ,t s ) (6) 



G(x,t <- x s ,t s ) 



( 1 \ 3/2 A 
\4irK At 



e "o 1 



(7) 



where At = t — t s and (Ax) 2 = (x — x s ) 2 . The unique 
argument of the Green function is actually the diffusion 
length A = AKoAt, because the energy dependence enters 
only in this combination. This is the characteristic length 
of an electron traveling during its propagation. Here we 
are going to focus on diffusion models for which the half- 
thickness Lh ~ 4 kpc is small compared to the radius of 
the disk Rh , such that in practice the radial boundary has 
negligible effect on propagation [33] . The Green function 
for which n(x, E) vanishes at z = may be expressed 
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as 



n— oo 

G(A, L) = —j=—^ J2 (-l)«exp 



— (2nL + ( — l) re -s s ) /A 



(8) 

As for the source term, in the case of production from 
DM annihilation with cross-section (av), it can be ex- 
pressed as 

( w\ ( \ J p{x) V dN{E) M 
q{x,E) = rj (av) i } — . (9) 



m DM J 

Here p(a;) is the DM profile, given in units of GeV/cm 3 , 
and rj = 1/4 or i] — 1/2 depending on the Dirac or Majo- 
rana nature of DM. The injection spectrum of electrons 
is given by dN/dE. 

With all these ingredients we can express the electron 
number density as 



n(x, E) = rj (av) 



( Po 



\m DM J b(E) 



1 



mDM dN 



(10) 

Here /(A) is the so-called halo function, which is defined 
has 



= / d^sGix^E^x^Es 

'DZ 



P{Xs) 

Pq 



(11) 



When electrons and positrons are created in the Galaxy, 
they are accelerated by the local magnetic field, thus pro- 
ducing synchrotron radiation, with a flux, per unit fre- 
quency v per solid angle arriving at position x, that is 
given by 



F v = — J dl J dE e P{v) n(x,E e ). 



(12) 



where integration is performed along the line of sight 
(los). Here P(v) is the synchrotron energy loss per unit 
frequency (in GeV/s/Hz), and together with the n(x, E) 
it must be integrated over all electron energies E e that 
produce the same synchrotron at frequency v. 

In order to express P(v) as a function of the elec- 
tron energy, we use the conservation of energy, P{y)dv = 
(dN e /dE)P(E)dE. Here P(E) is an energy-loss rate in 
units of (GeV/s), and dN e /dE is the number of elec- 
trons with an energy comprised between E and E + dE 
(not to be confused with dN/dE in Eq. Q). It can 
be demonstrated that the spectrum of an electron gy- 
rating in a magnetic field has its prominent peak at 
the energy corresponding to the critical frequency v c , 
dN e /dE ~5{E~E C ), where 



So 



<IE, /'(>.) mr.E, ) - ( ^Lp(E e ) n{x,E e ) 



with which finally 
1 



F„ = 



dl[ d ^P(E)n{x,E) 
4tt 7ios V dv 



(14) 



E=E C 



Deriving Eq. (13) with respect to frequency and using 



Eq. (10) and Eq. (14), gives the final expression: 



F„ = 1.21 x 10' 



;Jy (v 



sr [2 V3.1 x 10" 26 cm 3 /s 



(15) 



x Z lGcV/c 2 
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P& 



lGeV/c 2 /cm 3 



1/iG 
B 



'1GHz P(E C ,B) 1 f dl f mDM , dN , 
W / — / dE s —I{E,E s ) . 



v b{E Cl B) 4tt J kpc 



dE' 



where the synchrotron flux is expressed in Jansky (Jy), 
Uy = 10~ 26 W-m _2 -Hz _1 In this expression, the 
synchrotron loss rate is given by 



P(E,B) 



e 4 E 2 B 2 



(16) 



^-(uiv) (£) g « v -' 

The total energy loss assumed here is given by 

b(E,B)=P(E,B)(l + r lc/syn ) (17) 

where 

-2 



r IC/syn 



2 U ra d 



B 



3 B 2 /2 / i 



8 eV/cm I V 10 P G 



is the ratio between IC and synchrotron energy loss and 
t^rad is the total radiation density. 



B. Numerical approach 

In parallel to the semi-analytic expressions, we have 
determined the synchrotron fluxes using a fully numer- 
ical approach as implemented in the GALPRDP v54 code 
[HI [2D]. We ran GALPROP in 2D in galacto-centric cylin- 
drical coordinates (i?, z), solving the CR transport equa- 
tion on a grid and assuming cylindrical boundary condi- 
tions. Both re-acceleration and convection are included 
in GALPROP, though we set V c to zero in our calcu- 
lations for simplicity. The energy losses include also 
bremsstrahlung losses, based on up-to-date whole sky 
maps of the interstellar gas (for details on generation 



4 In the following, we will express the flux F v in terms of 
the brightness temperature of the radiation, T, where F v = 



E F —E r 



[cxp(hufkT) — 1] 1 is the usual black-body relation. 
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of gas maps see [22]). For IC losses GALPROP uses 2D+1 
map^Jof the interstellar radiation field, computed based 
on a model of the radiation emission of stellar popula- 
tions and further reprocessing in the Galactic dust [TP] . 

The spectrum and distribution of the synchrotron 
cmissivity as seen by an observer at the solar position 
depends on the form of the magnetic field, and the spec- 
trum and distribution of CR leptons and is computed as a 
function of [x , y, z, v). The emissivity is then integrated 
over the line-of-sight to get the synchrotron intensity. We 
chose the form of the magnetic field consistently with the 
semi-analytic approach, as detailed below. 



III. GENERAL ASTROPHYSICAL SETUP 
A. Astrophysical uncertainties 



Model 


Lh [kpc] 


K [cm 2 s" 1 ] 


S 


v a [km s 1 ] 


MIN 


1 


4.8 10 26 


0.85 





MED 


4 


3.4 10 27 


0.70 





MAX 


15 


2.3 10 28 


0.46 





la 


4 


6.6 10 28 


0.26 


34.2 


lb 


4 


6.6 10 28 


0.35 


42.7 


2a 


10 


1.2 10 29 


0.3 


39.2 


2b 


10 


1.05 10 29 


0.3 


39.2 


PD 


4 


3.4 10 28 


0.5 






TABLE I. Upper three raws: parameter sets derived using a 
semi-analytical approach, to lead to MIN/MED/MAX anti- 
proton fluxes at Earth from an exotic Galactic component 
[21] . Lower four rows: parameters from [23], consistent with 
[22] , derived using GALPROP code in a fit to CR data, and 
shown to reproduce the gamma-ray diffuse data well. Last 
row: plain diffusion model, shown to be consistent with the 
radio data at 22 MHz - 94 GHz frequencies, [24]. 



Studies coming from N-body simulations have led to 
popular expressions for the distribution of DM in the 
Galactic halo, like the Navarro-Frenk- White (NFW) den- 
sity profile [33] 



r 



PnfwM 



or the Einasto profile 



l>E\\(r) = /».,c.xp <j ((~ I ~ 1 



(19) 



(20) 



where r the radial distance from the centre of the DM 
halo. On the other hand, observations of galactic rotation 
curves as well some of the simulations which include also 
baryonic feedback on DM density (for a recent example 
see [37]) find DM density profiles which are more cored 
towards the inner regions of the Galaxy. One example of 
such profile is the isothermal profile [35] [35] 



Piso(r) 



(21) 



or modified Einasto profile (in this case the parameter a 
is smaller when compared to the parameters found in sim- 
ulations which contain only DM component, [3D] ). While 
the parameter a for the Einasto profile is fixed from a fit 
to the simulations, the values of parameters p s , a typical 
scale density, and r s , a typical scale radius for the Milky 
Way are determined from astrophysical observations (e.g. 
local stellar surface brightness, stellar rotational curves, 
total Milky Way mass within a given distance..., see e.g. 
[HUD]). As we will see in section 11IC the ratio of the 



Two spatial and the frequency dimension. 



synchrotron signal calculated with DM density of these 
three profiles gets smaller at higher latitudes, as at those 
distances (closer to the Solar position) DM density is bet- 
ter constrained. As the rotational curve measurements 
are poor at distances smaller than 2 kpc (or ~ 10°) from 
the Galacic Center and those regions of the Galaxy are 
baryon dominated, we choose the Region Of Interest 
(ROI) which spans |6| £ (10 ±3)° in Galactic latitudes, 
and | Z | < 3° in Galactic longitudes. 

Together with DM density profile, the CR propagation 
parameters pose one of the main uncertainties in predic- 
tion of the synchrotron signal. In [16] . the three models 
called MIN/MED/MAX [21] and featured in Tableware 
used to probe the uncertainty in CR propagation param- 
eters. Originally, those parameters were derived to pro- 
duce the maximal, median and minimal anti-proton flux 
from dark matter, while being compatible with the CR 
secondary to primary B/C ratio measurement [32 132] • 
Therefore, by construction, they do not necessarily cap- 
ture the uncertainty in the electron fluxes in the inner 
Galaxy, which is of interest here. In section [ill D| we will 
comment in more detail on the impact of a choice of a 
CR parameters, exploring additional sets consistent with 
the CR data and which were i) derived numerically and 
ii) shown to reproduce the observed whole sky gamma- 
ray or radio emission, therefore probing more directly the 
signals in the inner Galaxy, [2"2Tl24j . 

The Galactic magnetic field (GMF) is considered pos- 
sibly the most important ingredient when dealing with 
synchrotron radiation. In the diffuse interstellar medium 
it has a large-scale regular component as well as a small- 
scale random part, both having a strength of order micro- 
Gauss. The best available constraints in determining 
the large-scale GMF are Faraday rotation measures and 
polarized synchrotron radiation, see e.g. [4"3M46| . while 
random component is deduced mainly based on the syn- 
chrotron emission. 



G 



Several 3D models of the large scale magnetic field are 
implemented in GALPROP (e.g., [4TJ [48] ) . However, as our 
ROI is away from the Galactic disk and the small-scale 
random component is expected to dominate the total 
value of the magnetic field, we use a simple parametrisa- 
tion for a total magnetic field as customary in literature: 



B{p 1 z) cx exp — 



(22) 



The parameters R m and L m should in principle depend 
on the diffusion model assumed (or vice versa) , since the 
propagation in the Galactic medium is intimately related 
with the magnetic field. We follow [TB] in taking L m — 
6 ■ Lh and R m — 5 ■ Rh- It was shown in that work 
that an actual extent of L m and R m does not play a 
critical role (the difference between a constant magnetic 
field and the exponential form defined above is <30%), 
as long as the field extends into our ROI (i.e. L rn > 1 
kpc). The normalisation at Sun's position Bq is more or 
less well constrained, and we take the value of 6 /iG [4"5] . 
consistent with the measurements [43rl46| . However the 
value of the field in the inner Galaxy is cons ider ably less 
known and we rewrite the normalisation of ( 22 1 as [^] 



B a = B G [l + K G(R IG -r)] 



(23) 



where <3(Rig — r) is the unit step function as a function 
of r = \J p 2 + z 2 . With this change, while leaving Bq 
unchanged locally we allow for the magnetic field to have 
a higher effective normalization Bq (1 + K) in the inner 
Galaxy. We arbitrary set Rig = 2 kpc, in order to cover 
our ROI. We therefore fix the spatial dependence of the 
B field and explore the impact of overall normalisation 
Bq (1 + K) in the inner Galaxy, on synchrotron fluxes in 
section IIIIEI 

Contributions to the energy losses for electrons are as- 
sumed to come only from synchrotron and IC processes 
in our semi-analytical approach, as commented above. In 
principle, the radiation density i7 rac i (see Eq. (18)) has a 
spatial profile, which affects the synchrotron flux estima- 
tions. However, in the semi-analytical estimations used 
here we take t/ ra( j to be constant, which turns out to be 
a good approximation for the particular study performed 
in this work, as is justified in the next section. For the 
value of t/ ra d we take 8 eV cm -3 , which is approximately 
the value read off the ISRF maps used by GALPROP in our 
ROI.Q 



Concerning the frequency of observation we focus on 
the data taken at 45 MHz. In [16] it was shown that in a 
wide range of frequencies (22 to 1420 MHz) the change in 
the measured synchrotron flux is very small, while going 
to lower frequencies maximises a synchrotron signal of 
low mass WIMPs. The RMS temperature noise of this 
survey is 3500 K, which is generally subdominant with 
respect to the systematics errors involved in theoretical 
modeling. 



B. Cross-check of semi-analytical calculation of 
synchrotron fluxes 

In Fig. [T] we show a comparison of synchrotron spec- 
tra obtained with the semi-analytical (SA) and numeri- 
cal/GALPROP (NG) approaches. We see in general a good 
agreement between the two predictions. We observe that 
in the region of low latitudes, SA calculation gives rel- 
atively larger fluxes than the NG counterpart, while for 
the region of large latitudes, it is the opposite. That is 
expected, as in the SA approach we have used a constant 
value for the radiation density, while in the NG 
full space dependent profile for U la d is used. Typically 
this radiation density gets smaller as the distance from 
the galactic centre increases, so the synchrotron flux be- 
comes less suppressed, explaining why the GALPROP es- 
timations are typically larger for large angles and vice 
versa, see Fig. [T] 

However, as commented in section D, the data of syn- 
chrotron flux we used to constrain dark matter models is 
the one taken at a latitude of ~ 10°, and for this latitude, 
the semi-analytical and numerical approaches are in very 
good agreement. In other words, for this latitude the 
approximation of taking C/ ra d constant works very well. 

Finally, we have checked a good agreement with the 
results of 16] using their particular setup. The biggest 
difference between our works is that we assume a factor of 
~ 2 higher ISRF energy densities, inspired by the values 
implemented in the GALPROP code in our ROI. That in 
turn implies higher overall energy losses, and therefore 
lower synchrotron fluxes (by about the same factor, if 
IC losses are the dominant ones, see Eq. (24)), for the 
otherwise same parameter set. 



6 following [501 . 

7 Note however that the conditions in the inner Galaxy might 
be quite different from a simple CR propagation setup assumed 
here. In particular, observations of the bubbles like structures 
Centerd at the Galactic centre and extending to 50 deg in lat- 
itudes in gamma rays |51l I52| and WMAP haze at microwave 
frequencies, observed by WMAP |53j and Planck |54| . witness of 
possibly more complicated configuration of the magnetic fields 



and CR propagation parameters in that region. While bubble- 
like structures appear sub-dominant with respect to the standard 
components of the diffuse emission in our ROI and electron ener- 
gies of interest, we caution that before their origin is understood, 
the actual structure of magnetic fields or the CR propagation 
conditions, cannot be reliably modeled. 
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FIG. 1. Comparison between semi-analytical (S-A.) and numeri- 
cal/Galprop (N-G) approaches, for different annihilation channels, 
and DM mass of m jjm = 10 GeV. A frequency of 45 MHz was used, 
assuming a magnetic field normalised to Bqc = 10 l*G (which cor- 
responds to K=0, i.e. no enhancement of the magnetic field in our 
ROI). 



Synchrotron signal for different choices of DM 
density profile 



In the remainder of the text we will focus on two DM 



profiles, NFW, eq. (|19j) and isothermal (ISO) Eq. ([21} , 
using the following values of parameters, consistent with 
observations: p s — 0.31 GeV cm -3 , r s = 21 kpc, for an 
NFW profile, and p s = 1.53 GeV cm -3 and r s = 5 kpc 
for Isothermal profil^J [16] . However, in this section we 
also show the prediction for the synchrotron signal in 
the case of a modified Einasto profile. In particular, this 
profile is modified to have a shallower inner slope than the 
usual Einasto profile found in DM-only simulations and 
it describes better results of simulations which include 
baryonic feedback (parameters we use are a = 0.11, r s = 
35.24 kpc, p s = 0.041 GeV cm~ 3 [40 ). As this profile is 
'bulkier' at distances 1 kpc from the GC, the DM signals 
are generally higher than those of NFW in that region. 

In the parameter sets we chose, the local value of DM 
density is set to p G — 0.43 GeVcm~ 3 . One should also 
keep in mind that the overall normalisation of DM distri- 
bution pq is uncertain, being in the (0.43±0.113±0.096) 
GeV cm -3 rang^] [55] . Therefore, in addition to the dif- 
ferences in the signal caused by the DM profile shape, 
and shown in Figure [2j synchrotron signals scale with 



8 Note that we make a conservative choice by choosing a rather 
extended core. Smaller values of r s would result in fluxes more 
similar to those obtained with the NFW profile. 

9 However, depending on the analysis, the uncertainty window on 
this value can vary in the 0.2 - 0.8 GeV cm -3 range. 
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FIG. 2. Comparison between synchrotron signals for a DM mass of 
10 GeV annihilating to muons (top) and b quarks (bottom figure) 
for three DM density profiles: modified Einasto, NFW and ISO 
thermal profile. pQ = 0.43 GeVcm~ 3 is assumed for this plot (see 
text for the remaining parameters) and propagation of electrons is 
done using a CR propagation setup as shown in Table [I] 



D. Synchrotron signal for different choices of CR 
parameters 



As discussed above, MED /MIN /MAX sets of CR prop- 
agation parameters were derived using a semi-analytical 
description of CR propagation, and a fit to B/C measure- 
ment, with a requirement to produce minimal, medium, 
and maximal DM generated anti-proton fluxes, at a Solar 
position. These models were recently reanalysed in [56] . 
where whole sky radio data together with B/C measure- 
ments where considered, testing therefore impact of these 
parameters on electron population, and its synchrotron 
emission in various regions in the Galaxy. This work con- 
cludes that both MIN and MAX models are disfavored 
by radio data towards the galactic anti-centre. Demand- 
ing consistency with B/C nuclei data, it is found that 
small halo sizes ~ 1 kpc are essentially excluded (see 
also [21]), but also large values Lh > 15 kpc show some 



tension with radio data, arguing that MIN/MAX sets 
of parameters present somewhat extreme choices when 
compared to observations complementary to CR nuclei 
data. 

In parallel to the above mentioned analysis based 
on semi-analytical approaches, we also use results of 
three analysis based on a numerical calculation with the 
GALPROP code [22H24] . In [22] the GALPROP code is used 
to derive a set of CR parameters which provide a good 
fit to the CR data, and at the same time reproduce the 
gamma-ray data well. In a parallel work, [23 , a full 
Bayesian analysis of a fit of CR models is made using 
GALPROP to confront CR data and to derive the best fit 
model and its scatter. It is shown that gamma-ray pre- 
dictions of such model are consistent with the Fermi LAT 
gamma-ray measurement. We will use the best fit model 
from |23j , and few models within a 2 sigma scatter of the 
Bayesian analysis, see Table [I] We show a comparison 
of synchrotron signal calculated with these choices of CR 
propagation parameters in Fig[3] The GALPROP based pa- 
rameter sets in general have lower values of parameters 8 
and consequently higher values Kq are obtained in a fit 
to B/C data. That in turn means that diffusion is higher 
in this set of models, i.e. electrons diffuse out of our ROI, 
and contribute instead more at higher latitudes, as seen 
in Fig. 

A study in [53] has shown, however that some models 
with re-acceleration (v a ^ 0, similar to the case la above) 
are in tension with the large set of radio data. The main 
reason for this is that production of CR secondary elec- 
trons is enhanced in this case and produces too strong 
radio emission at low frequencies. We therefore consider 
also a plain diffusion model (PD) often used in literature 
and shown to be consistent with the radio data in |24j . 

In the remainder of the text we will adopt the MED 
model together with the NFW profile and p = 0.43 GeV 
cm" 3 as an moderately optimistic scenario, and the 'la', 
GALPROP based CR propagation model, together with the 
Isothermal profile and p© = 0.3 GeV cm~ 3 as a moder- 
ately conservative setup. 



E. Synchrotron signal for different magnetic field 
choices 



In this section we study how DM limits change de- 
pending on the assumptions of the overall normalisation 
of the magnetic field in our ROI. It has already been 
noticed (see e.g. [TS]) that for a fixed electron injection 
spectrum there exists an optimum value for the magnetic 
field which maximises the synchrotron flux at a given syn- 
chrotron frequency. However in this section we want to 
understand this fact in more detail. 

As can be seen in Eq. [15] the magnetic field influ- 
ences the flux through the energy losses, and through 
the electron energy. In principle, the electron spectrum 
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FIG. 3. Comparison between synchrotron signals for a DM mass of 
10 GeV annihilating to muons (top) and b quarks (bottom figure) 
for different CR propagation setups, detailed in Table]!] NFW DM 
profile is assumed here. 



dN(E s )/dE affects the maximisation of the flux with 
respect to the magnetic field B. A way to see this is 
the following: in the assumption of considering a one- 
to-one correspondence between the emitted synchrotron 
frequency v and the electron energy E c , (see Eq. 131, for 



a given value of the magnetic field B, the correspondent 
E c could be disfavored (or not available) by the elec- 
tron spectrum produced by the DM. However, as in this 
study we are going to work with synchrotron data with 
frequency of 45 MHz, and magnetic fields > 6 /xG, the 
electron energies producing those desired frequencies are 
always smaller than 1 GeV, which can perfectly be pro- 
duced by our DM candidates in the range [1-200] GeV. 
So in practice, for our analysis, the only dependence on 
B that matters are encoded in the energy losses. Indeed, 
given a specific frequency v, and a given annihilation 
channel, it can be expressed in the following form: 



F(v, B) oc 



B 2 



1 



(24) 



u + B 2 ) y/B 

where a represents here the rest of energy-losses, here 
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assumed to be only IC. Note that since both synchrotron 
and IC losses scale with energy as E 2 , the energy depen- 
dence cancels in this particular analysis. From (24) one 



observes two extremal cases: one in which synchrotron 
loss is negligible with respect to the rest of energy losses, 
for which the flux scales with B as F ~ B 3 ' 2 , thus in- 
creasing as B increases, and the other, in which syn- 
chrotron is actually the dominant energy loss, after which 
the flux scales as F ~ 1/vB, thus decreasing as B in- 
creases. In other words, there will be an intermediate 
value of the magnetic field for which synchrotron becomes 
the dominant energy loss, and this value is actually the 
one maximising the flux. 

Figure [4] shows the shape of synchrotron flux as a func- 
tion of the value of magnetic field at GC. Taking into 
account only IC (apart from synchrotron of course), one 
can have an idea about the maximum of the flux already 



by direct differentiation of ( 24 ) , assuming the values of a 



correspondent to this case. The value of B for which the 

For U lad = 8 



flux is maximal scales as Bq^ 



eV/cm 3 , B^ x ~ 26 juG. 

In the following we will show DM limits using two 
cases for magnetic field configurations. The first one, 
a standard Bqq — 10 /xG value in our region of in- 
terest is assumed, and the second case will be the nor- 
malisation which maximises the synchrotron signal, i.e. 
B GC ~26fiG. 



Bqc (hQ) 



FIG. 4. Flux predicted at 10 deg off the GC, as a function of 
the value of magnetic field at the GC (roughly, at our ROI, we 
have -Broi ~ 0.5-Bqc)- Lines represent the results in the Semi- 
Analytical approach. We assume a DM annihilating directly to 
electrons, using the NFW profile and the MED diffusion model 



IV. OTHER INDIRECT CONSTRAINTS 

In the following we will compare the limits we derive 
from considering radio data to i) those derived by the 
Fermi LAT collaboration from the non-observation of 
dwarf spheroidal Galaxies (dSphs) in gamma-rays |17j 
and ii) constraints based on the measurement of the 
CMB anisotropy spectra [IS]. These limits are among 
the strongest to date and uncertainties involved in their 
calculation are complementary to the ones derived in this 
work. 

The main strength of considering dSphs is that they are 
DM-dominated and lack active astrophysical production 
of gamma-rays. We will refer to an analysis of the Fermi 
LAT data which, for the first time, combines multiple 
Milky Way satellite galaxies in a single joint likelihood 
fit and includes the effects of uncertainties in J factors 
(line of sight integrals of DM density squared, which de- 
fine the DM annihilation rates). The improvements of 
the limits derived in such an analysis over 10 years of 
Fermi LAT mission was estimated in [ST]. In the low- 
energy regime, the sensitivity increases as roughly the 
square-root of the integration time. However, in the high- 
energy (limited background regime) the LAT sensitivity 
increases more linearly with integration time. Thus, 10 
years of data could provide a factor of V5 to 5 increase 
in sensitivity. Additionally, new optical surveys (such 
as Pan-STARil^Jand the Dark Energy SurvejQ could 
provide a factor of 3 increase in the number of detected 
dSphs corresponding to an overall increased constraining 
power y/lE to 15. As we focus mainly on DM masses 
< 100 GeV, we will use y/15 to illustrate potential im- 
provement in our plots. 

In order to translate these limits in terms of democratic 
couplings to leptons and hadrons (shown in Figures [6] 
and 7| we use i) the limits to r channel worsened by a 
1 /3 branching to this channel; dSphs constraints on (x are 
significantly weaker (this likely holds true also for the e 
channel, however these limits were not published) and ii) 
limits on b quark states, as the gamma-ray spectra are 
very similar for all quark channels. 

CMB constraints arise from redshifts in the range 
100 < z < 1000 Ql [58H5T] . The physical effect of energy 
injection of (exotic) particles around the recombination 
epoch is that it results in an increased amount of free elec- 
trons, which survive to lower redshifts and increases the 
width of last scattering surface, consequently suppressing 
the amplitude of some of the oscillation peaks in the tem- 
perature and polarisation CMB power spectra. As con- 
straints come from high redshifts well before the forma- 
tion of any sizeable gravitationally bound structure, this 
set of constraints does not depend on highly uncertain 



10 http:/ /pan-starrs. ifa.hawaii.edu/public/ 

1 1 http: / / www.darkcnergysurvey.org/ 
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parameters related to structure formation. Detailed con- 
straints have been recently derived in |18j . based on the 
WMAP (7-year) [52] and Atacama Cosmology Telescope 
2008 data [63 . The constraints are somewhat sensitive 
to the dominant DM annihilation channel: annihilation 
modes for which a portion of the energy is carried away 
by neutrinos or stored in protons have a lesser impact on 
the CMB; on the contrary the annihilation mode which 
produces directly e + e~ is the most effective one. The lim- 
its on hadronic channels were not derived in this work. 
We therefore show CMB limits only for the democratic 
coupling to leptonic states, in which case we combine the 
published limits to /i and e channels, rcnormalised by a 
branching factor of 1/3. We also show the projected sen- 
sitivity of the PLANCK telescope in the near future, as 
calculated in [T51. 



CONSTRAINTS ON DARK MATTER 
MODELS 



data produced at LEP [5] and put lower bounds on the 
effective coupling A e . This effective approach can be un- 
derstood as a limit of microscopic models with heavy 
mediator (Higgs boson or extra t/(l) Z' for instance), 
discussed in sections IV Bl and IV CI 

A similar study has been made [5J [7J using results from 
the CDF collaboration with Tevatron [57], and more re- 
cently [9\ from ATLAS study [68], and by the CMS Col- 
laboration itself [TU]. The analogous types of operators 
defined in Eq. [25] are used as source of processes like 
(qq — > jxx)i where q are SM quarks and j stands for 
a single jet. A mono-jet final state is then used to put 
constraints on A q and thus on (in this case) hadronic 
annihilation cross-section, in the very same way as in [5] . 



scalar, synch. (lep).26iiG ■ 
scalar, synch. (lep).10uG - 
vector, synch. (had). 26uG ■ 
vector, synch. (had). 10uG - 
scalar, LEP • 
vector, CMS(had) ■ 



A. Effective operators 



Mono-events at colliders 



CD 




The very same interactions responsible for DM anni- 
hilation in the galactic medium that ultimately can pro- 
duce synchrotron fluxes, are the ones by which a DM 
signal could be produced and measured in collider exper- 
iments. At LEP, for example, when a pair of electron- 
positron collides at high energy a plethora of final states 
are produced. Among these states, the production of 
a DM candidate \ can be characterised by missing en- 
ergy and the emission of a single photon [5H |5S] . These 
data can then be used to constrain theoretical models. 
We first choose to use the powerful machinery of effec- 
tive field theory (EFT) to capture features from a broad 
class of WIMP models. In particular we assume that at 
the relevant energies, the only available degrees of free- 
dom are the dark matter particle itself and the Standard 
Model. Provided this is true for all energies of interest, 
the EFT provides a common language which allows one 
to compare the constraints from different types of exper- 
iments [55] . Introducing new physics scales A, one can 
define effective couplings between {e.g. fermionic) DM 
to SM fermions / by 

0{ = —mxXf^f), 0{ = ~( X f)(h), (25) 

where T; =(1;7^; 757^) for scalar (OcJ), vector (O v ) or 
axial (0^) operators respectively. Performing a Monte 
Carlo simulation studying the signal (e + e~ — » 7xx) phis 
background (e.g. e + e~ — > jDiy) processes, one can com- 
pare the theoretical background predictions with the real 
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FIG. 5. Lower bounds on the effective scales Ai q from some col- 
lider studies, as well as from synchrotron data for a certain choice of 
parameters (see text for details). "CMS(had)" means CMS bounds 
to quark operators. Solid lines are bounds from synchrotron radia- 
tion data, while dashed lines are bounds from LHC (CMS) or LEP. 
Red lines correspond to quark operators while black lines corre- 
spond to leptonic ones. 

As an illustration, we show in Fig. [5] the limits we 
obtained on the leptonic and hadronic effective scales, 
A; and A q , obtained by LEP and CMS as function of 
the DM mass m x in the case of scalar /vectorial interac- 
tion with universal leptonic/hadronic coupling. We no- 
tice that the bounds on A^ = ; q become weaker as DM 
mass increases because of the phase space reduction, up 
to some point in which it is kincmatically impossible - 
given the energies of the experiment- to produce the pair 
of DM particles (the threshold is of course higher in LHC 
experiments than at LEP). The lower bounds on A;(q) 
is around 300 GeV (600 GeV) and relatively indepen- 
dent of the DM mass as one measures missing energy. 
It becomes then interesting to compare this limit to the 
one derived from synchrotron emission. These bounds 
obtained by the LEP and CMS can also be converted 
directly into upper bounds on the leptonic/hadronic an- 
nihilation cross-section {ov) ee / qq from the inverse pro- 
cesses (xx — > ee/qq), generated by the operators defined 
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in Eq. 25 (see Appendix |A| for details). 



Synchrotron vs Collider and complementary Indirect 
Detection bounds 



In order to model the synchrotron signal of a standard 
astrophysical origin, a HASLAM 408 MHz map jSHj can 
be used as a template. It is a full sky radio map with 
the best angular resolution and sensitivity. Indeed, a 408 
MHz frequency is well suited to gauge the contribution 
from the harder population of galactic electrons, while 
leaving lower frequency maps sensitive to test a possible 
contribution from softer electrons originating in annihi- 
lations of light dark matter candidates. 

In |16| it has been shown that by extrapolating the 
HASLAM data down to 45 MHz, to model astrophysical 
synchrotron emission, one is left with ~ 30% residuals 
when compared to the actual 45 MHz data. Those resid- 
uals do not have the proper morphology of a DM signal 
and therefore no DM detection could be claimed. Based 
on this analysis it seems reasonable to assume that the 
current uncertainty (i.e. systematics) in modelling the 
astrophysical emission is at a level of ~ 30%. However, 
in our analysis, we do not attempt to model specific as- 
trophysical signals. Instead, we decide to apply a con- 
servative 95% CL DM limits without assuming any con- 
tribution from astrophysical background, and which we 
label no bckgd, in the figures and also show how this lim- 
its change as a function of the systematic uncertainties 
on the astrophysical model, expressed as a background 
uncertainty in % of the data. 

Our procedure is to compute the synchrotron flux com- 
ing from DM (Fdm) in our ROI, and then constrain the 
result with the available data corresponding to that re- 
gion, F bs, by requiring that _F DM < F ohs + 2a, where 
la- is the uncertainty considered in each case (see Fig- 
ures [8] and [9]). In the case of no bckgd, la is the rms 
temperature noise, taken to be 3500 K. 

For illustration, in Figure [5] we have included the 
bounds on Aj obtained from synchrotron data at 45 
MHz with a 26 and 10 fj,G magnetic field in our ROI, 
with the hypothesis that the dark matter signal lies 
within 5% of uncertainty of the background, and using 
the NFW+MED set-up. We clearly see that the syn- 
chrotron radiation, with the present data, can already 
give stronger limit on the effective scales, independently 
on the value of the magnetic field, and is complementary 
to the bounds from accelerator searches. Indeed, whereas 
there are no "threshold effect" for synchrotron radiation 
at large DM mass, one notice that for low masses, in 
the hadronic channel the bounds becomes weaker. We 
comment on this behaviour below. 

We also translated these limits on Aj into limits on anni- 
hilation cross sections. We present the result in Figures [6] 
and [7] for different values of the magnetic fields, and na- 
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FIG. 6. Synchrotron bounds on {uv), assuming DM couples 
democratically to all charged-leptons (left panel) or to all kinemati- 
cally available quarks (right panel). In both cases, comparison with 
bounds coming from colliders are shown explicitly. The magnetic 
field normalisation has been set to 10 ^G. Astrophysical setup: 
MED diffusion model, NFW profile with p e = 0.43 GeV/cm 3 . 
Purple lines represent present (solid line) and 10-years projection 
(dot-dashed line) limits coming from Dwarf galaxies. Brown line 
in the top panel represents current DM limits based on the mea- 
surement of CMB anisotropies (solid) and near future reach based 
on the Planck data (dashed line). 



tures of the DM couplings, in comparison with limits 
coming from LEP, CMS, observation of dwarf galaxies by 
Fermi LAT and the measurement of CMB anisotropies 
by the WMAP and ACT. Comparing the 2 figures, we 
first notice that the results are not extremely sensitive to 
the value of the magnetic field (10 or 26 /iG) and that 
synchrotron bounds largely compete with bounds from 
other type of experiments. Secondly, the hadronic chan- 
nel gives the weaker bounds on (<tv), especially at low 
masses. For low masses, the annihilation into hadronic 
states gives an electron spectrum much softer, especially 
after the kinematic closure of the bb channel, and ra- 
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FIG. 7. Same as in figure fill, with the same astrophysical setup, 
but now the magnetic field normalisation has been set to 26 /J.G. 



diates much less synchrotron emission than in the case 
of a pure leptonic channel: a lot of soft electron pro- 
duced by later decays of final state hadrons will radi- 
ate at frequencies lower than 45 MHz and will not be 
bounded by current radio data. At the high mass end, the 
bounds coming from indirect detections searches (based 
on the Galactic synchrotron emission, dwarf Galaxies or 
CMB) get weaker because the flux of products of DM 
self-annihilation is proportional to 1/mt. We should 
also emphasise that even the very conservative result, 
where we supposed that all the data are generated by 
DM synchrotron, gives already limits competitive with 
LEP bounds, whereas the 5% uncertainty bounds could 
give the best limit obtained by an indirect detection ex- 
periment. 

Concerning the accelerator constraints, from the study of 
the (<jv) predictions for different kind of effective interac- 
tions (cf. Appendix [A]) , one realizes by direct inspection 
of those expressions that normally the vector operator 



gives larger values, with respect to the scalar and axial 
case. This is because its Lorentz structure allows for (o~v) 
to develop a term independent on v 2 and proportional 
to m^., and thus not suppressed. The consequence of 
this is that in general, vector interaction produce weaker 
bounds. 

In the comparison with accelerator constraints, one re- 
alises that in the case of leptonic channel, the synchrotron 
bounds in the conservative case of assuming no back- 
ground, (orange band) is in general as competitive as the 
bound coming from LEP studies, assuming a vector ef- 
fective interaction. It starts to be more competitive for 
masses larger than tens of GeV, and beyond 100 GeV, 
space of parameter space beyond the reach of the accel- 
erator. In the case of CMS bound from a vector interac- 
tion, it is all the way slightly better (for light masses) or 
similar (for heavier masses) than synchrotron constraints, 
obviously up to the TeV range where LHC starts loosing 
sensitivity. The LEP bound coming from a scalar effec- 
tive operator is however one order of magnitude stronger 
for light masses of ~ 10 GeV. Furthermore, to have an 
idea of how well we would need to understand the back- 
ground in order to render synchrotron more competitive 
than LEP or CMS, we include in Figure [6] a bound as- 
suming that background is known within 5%, at 95% CL. 
We observe that for these already small uncertainties the 
synchrotron searches start to be more competitive (in 
general) than collider bounds, independently of the DM 
mass (in the case of LEP). 

We also show in Figure [6] the bounds on hadronic chan- 
nels, now compared to collider bounds on vectorial and 
axial interactions. Again, for the conservative case of no 
background, synchrotron bounds start to be more com- 
petitive than collider bounds for tudm °f about tens of 
GeV, in the case of vectorial interaction, for ATLAS as 
well as for CMS. In the case of axial coupling, collider 
bounds are always stronger. In the optimistic case of al- 
lowing at most 5% of uncertainties in the background, 
the bounds from synchrotron can exclude those from col- 
liders already at 10 GeV in the case of vectorial interac- 
tion. The difference in the behaviour for masses lighter 
than 5 GeV, with respect to the leptonic case, is due to 
quark masses thresholds. In the hadronic case, when the 
bb channel opens, it gives an important contribution to 
the flux, causing the maximum allowed av value to de- 
crease. Since there no more threshold after the bottom 
mass (the top is heavier than the ranges considered here) , 
the flux decreases smoothly when tudm increases, as in 
the the case of leptonic channels (for which all thresholds 
are below ~ 2 GeV). 

We run a similar analysis in Figure [7J with the "opti- 
mum" choice of magnetic field, Bqq — 26 //G, studied in 
the previous section. We see in this case an improvement 
of synchrotron bounds with respect to collider bounds. If 
the DM interaction is vectorial, the synchrotron can ex- 
clude the collider bounds for all the mass range, for the 
case of leptonic channels, even in the conservative back- 
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ground assumption. For hadronic channel the exclusion 
is already effective for masses beyond 10 GeV. As be- 
fore, depending on our knowledge of the background, syn- 
chrotron could be able to exclude current collider bounds 
on axial interaction. 

To summarise, our analysis shows that synchrotron 
bounds, beside being complementary, could be poten- 
tially stronger than those obtained from collider searches, 
and that extensive studies of the astrophysical back- 
ground at these frequencies is well motivated. To this 
purpose consistent progress is expected to be achieved in 
the next years with the new high quality data coming 
from the PLANCK mission^] and from low frequency 
arrays like LOFAE0 which will survey low frequen- 
cies 10-240 MHz and is the pathfinder for the future 
SKA facility 

It is also worth noticing that the anal- 
ysis of the Galactic diffuse emission of the Fermi LAT 
will bring further insights on CR production and prop- 
agation processes, additionally constraining the Galactic 
synchrotron emission (for a current status see, [22" 1 124 ) ). 
We leave such dedicated studies for future works. 

The impact of the control of the background is illustrated 
in Figs. [8] and [9] where we show the level of precision 
that would be required to be able to exclude/probe the 
effective models for different DM masses. The dots on 
top of the iso-mass lines represent collider bounds for 
the corresponding mass. The lines with no dots corre- 
spond to cases for which the synchrotron limit is (for the 
whole range of uncertainties) stronger or weaker than the 
respective collider bounds. So, for each iso-mass, the re- 
gions on the left-hand-side of the dots correspond to sit- 
uations for which the synchrotron bounds are stronger, 
and on the right-hand-sides, to stronger bounds from col- 
liders. For instance, the case of a 20 GeV DM candi- 
date with scalar leptonic interactions is more constrained 
by radio synchrotron data, provided uncertainties on the 
background are smaller than 5%, than by the LEP data 
(Figure [8 - top). For vector-like interactions, the syn- 
chrotron bounds were always stronger, no matter what 
the uncertainties considered, for DM masses above 5 
GeV and assuming Bqq = 10 ^G (see Figure [8]). For 
Bqq = 26 fiG, (Figure [9]- top), the synchrotron bounds 
are always better than collider bounds, in the case of vec- 
tor interaction. The conclusion is that, for uncertainties 
on the background of 5%, or less, the synchrotron bounds 
are stronger than the LEP bounds, independently of the 
DM mass and of the nature of the effective DM coupling 
to SM leptons. 

Looking at hadronic channels, we see that for uncertain- 
ties of 5 — 10% (Figure [8 - bottom), the synchrotron 
bounds are stronger than those coming from LHC, if the 
interaction is axial-like and for Bqc = 10 ^G. The case 



12 http:/ /www.rssd.esa.int/index.php?project=Planck 

13 http://www.lofar.org/ 

14 http://www.skatelescope.org/ 
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FIG. 8. Bounds on (crv) coming from DM synchrotron fluxes 
as function of background uncertainties assumptions [in % of data] 
at 95%CL. Different choices for DM mass are plotted. Dots rep- 
resent the bounds coming from LEP (left-panel), specifically the 
case of a democratic coupling with leptons (as the synchrotron 
bounds) and assuming a scalar interaction; or LHC-ATLAS (right- 
panel), for democratic coupling with quarks, for both vector and 
axial interactions. The magnetic field normalisation has been set 
to 10 fiG. Astrophysical setup: MED diffusion model, NFW profile 
with p e = 0.43 GeV/cm 3 . 



of Bqc = 26 fiG (Figure [9] - bottom) is of course even 
better, and for m x ~ 60 GeV, even with background un- 
certainties of 20%, the synchrotron can rule out the col- 
lider bounds. If the interaction is vector-like, the bounds 
from colliders arc weaker, as wc discussed above, and for 
masses of 20 GeV and beyond, synchrotron bounds are 
always stronger. 
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dresses the abundance of WIMPs requires two fundamen- 
tal ingredients: 

• A candidate. It should be a massive and weakly in- 
teracting particle, stabilised by a (discrete or con- 
tinuous, possibly gauged) symmetry. 

• A mediator. This component can be an existing 
particle or a companion of the dark matter can- 
didate naturally present by construction. For in- 
stance it could be a particle (the Higgs, or a new 
scalar) or a new gauge boson (Z'). 

In this and the following section, we study the two sim- 
plest extensions of the SM with DM particles, the so- 
called Higgs-portal and Z', or kinetic mixing, portal. In 
particular we determine their radio synchrotron emission, 
and combined with the WMAP constraint, we analyse 
the parameter space that is experimentally allowed. 

The perhaps simplest extension of the SM consists in the 
addition of a real singlet scalar field, that couples to the 
SM through the Higgs field. In this case, the mediator is 
the Higgs boson itself, and the model is usually called the 
"Higgs portal" . Although it is logically possible to gener- 
alise this scenario to more than one singlet, the simplest 
case of one singlet already provides a very useful frame- 
work to analyse the generic implications of an augmented 
scalar sector. The most general renormalisable potential 
involving the SM Higgs doublet H and the singlet S is 



S 2 H^H 



H^HS 



-S* 



(26) 



FIG. 9. The same as in Figure [5] with the same astrophysical 
setup, but with the magnetic field normalisation set to 26 fiG. 



B. Higgs portal 

The effective operator approach is quite powerful, as it 
encompasses many possible underlying microscopic the- 
ories with a minimum set of free parameters, but it has 
its limitations. This is in particular the case when other 
degrees of freedom become relevant at the energies con- 
sidered (be in the early universe or at colliders). Also, a 
microscopic theory in general predicts some specific re- 
lations between the various effective couplings to which 
it reduces at low energies. Finally, a specific theory may 
have implications or predictions, that are hidden in an 
effective operator approach. For all these reasons, it is 
of interest to consider more fundamental models of dark 
matter. Hence the purpose of this section, which is to 
take our analysis of constraints from synchrotron radia- 
tion to a more microscopic level. 

Generally speaking, any fundamental theory that ad- 



For the case of interest, the S particle must be stable 
to be a dark matter candidate. This is achieved by im- 
posing a Z2 symmetry on the model, S — » —S, with H 
unchanged, thereby eliminating the K\ and K3 terms. We 
also require that the true vacuum of the theory satisfies 
(S) = 0, thereby precluding mixing of S and the SM 
Higgs boson (and the existence of cosmologically prob- 
lematic domain walls). In this case, the mass of the S is 
simply 



A, 



l s — Ms 



(27) 



with v — 246 GeV is the vev of the Higgs field, and the 
H-SS coupling is 



C 



HSS 



(28) 



where Mw is the W-boson mass and g the SU(2) gauge 
coupling. 

Different aspects of scalar singlet extension of the SM has 
already been studied in [70 whereas a preliminary anal- 
ysis of its dark matter consequences can be found in [71] . 
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Some authors considered the possibility of explain the 
DAMA and/or COGENT excesses [72], whereas several 
authors looked at the consequences of earlier XENON100 
data [72]. Other authors looked at the consequences on 
the invisible Higgs width at LHC [73] or the restriction 
of the parameter space due to an hypothetical 125 GeV 
Higgs signal [71] . Other studies probed the model by in- 
direct searches |75H77) of 7— ray or positrons, but there 
is not yet and analysis of constraints from synchrotron 
radiation. 

We show in Figure [To] (top) the synchrotron flux in our 
ROI expected for the Higgs-portal model for two different 
representative values of the magnetic fields (10 /iG and 
26 /iG) as a function of the dark matter mass, for a Higgs 
of 125 GeV and candidates that respect the WMAP con- 
straint (at 5a) . We also represented the limits we obtain 
based on the 45 MHz radio data in the two cases dis- 
cussed above: when we supposed that all the signal is 
due to synchrotron radiation (labelled "no backgr", ex- 
tremely conservative case A), and when we allowed the 
signal to lie within the 5% of uncertainties due to mod- 
elling of an astrophysical background of standard origin 
(labelled "5% background, less conservative case B). We 
do this for two typical astrophysical setups, as studied 
above in the paper: a NFW dark matter profile, with 
MED diffusion model and a local normalisation of a DM 
density of p & = 0.43 GeVcm -3 , and on the other hand, 
an Isothermal profile, with diffusion model " la" FM (see 
Table I) and a local normalisation of a DM density of 
p Q = 0.3 GeVcixr 3 . 

For the (NFW + MED) setup, we observe that, in the 
more conservative case, dark matter masses below ~ 20 
GeV are already close to be excluded for small values of 
the magnetic fields (B = 10 ^G), and even DM masses 
below 30 GeV are excluded for B = 26 p,G. This result is 
better that the one can obtain based on the Fermi LAT 
gamma-ray and PAMELA positron data [75j [76] , which 
is one of the important results of our analysis. Indeed, 
dwarf galaxies study exclude only m s < 10 GeV [77] . 
Also, with the same assumption on the CR propagation 
model (MED), one can see that the current synchrotron 
constraints already reach a higher level of precision than 
those based on the PAMELA positron data (see Fig. 2 
of [76]) and gives even better result than the projected 
AMS-02 sensitivity (see Fig. 3 of [76]) if the DM signal 
lies in the 10% uncertainty of the background. The ex- 
clusion by synchrotron emission is also similar to the one 
sets by the Fermi LAT data on galactic gamma-rays [75] . 

In case B, allowing the DM signal to stay within 5% 
uncertainty, we observe that mainly all the parameter 



15 Those models, where Alfven velocity can not be neglected - at 
least for our energies of interest-, can not be studied directly by 
the Semi- Analytical approach followed here. Instead, we use the 
full numerical analysis to estimate suppression factors in every 
annihilation channel, and from there we derive the new bounds. 



space of the model is excluded for masses lighter than 
100 GeV, except in two narrow regions, namely near 
the Higgs pole (2m s ~ m H ), and just above the W ± 
threshold {ms — 80 GeV). Indeed, in a region around 
the pole, the enhancement of the cross section due to the 
Breit-Wigner form of the amplitude implies that a very 
low value of Chss is required to match the WMAP con- 
straint, (av) ~ 3 x 10~ 26 cm 3 s _1 , as the amplitude is 
proportional to 



M cx 



R 2 



Chss 
2 ^iT 



in 



H 



with Ecm being the CM energy of the annihilating S 
pair. However, the synchrotron radiation produced in 
the Galaxy comes from S essentially at rest, whereas its 
kinetic energy at freeze-out was about T ~ m/25. Con- 
sequently, for annihilation in the Galaxy, the points re- 
specting WMAP are away from the pole and the enhance- 
ment is not sufficient to counterbalance the small value of 
Chss'- the annihilation cross section is (c^^ookms- 1 — 
10 _29 cm 3 s _1 . The synchrotron flux is thus largely re- 
duced in this region of the parameter space, explaining 
the first dip seen in Fig. [10] 

The second dip seen in Fig. [lO] is also due for kinetic 
reasons. Indeed, as soon as the W ± channel is open, the 
W final state is by far the dominant one due to the 
large gauge coupling of the W's to the Higgs. Thus, for 
singlet candidate with a mass slightly below Mw±, its 
kinetic energy at the freeze-out temperature allows its 
annihilation into W ± , even for small values oiCnss, a s 
the process is dominated by the gauge interaction and 
not the Yukawa ones (the bb final state is the dominant 
one for m s < Mw±)- However, nowadays, such slightly 
lighter singlet cannot annihilate into W ± , as its kinetic 
energy has dropped below the threshold. Consequently, 
its main decay in the Galaxy is into bb pairs, but with a 
small value of Chss, again giving a very low synchrotron 
emission. These two narrow regions are thus very difficult 
to observe with this type of signal. We should perhaps 
emphasize that the arguments we gave above are valid for 
any kind of, indirect detection prospect (7 or antimatter) 
and is not specific to synchrotron. 

We also show in Fig. 
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(bottom) the precision on 
the background that is required to be able to ex- 
clude/discover a synchrotron signal for a singlet DM 
through the Higgs portal. We observe that, except in 
the two narrow pole regions discussed above (where an 
unreasonable precision would be needed to measure any 
type of fluxes), if 5-10% (at 95% CL) of the observed 
flux is due to synchrotron radiation, the model is already 
excluded for Bqq = 26 pG. For Bqc = 10 ^G one needs 
to suppose that < 5% of the measurement is due to the 
annihilation of the singlet to exclude the model. This 
result is then essentially independent of the precise value 
of B GC - 

Of course, these conclusions are weaker for the second 
astrophysical setup considered here, where electrons dif- 
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fuse more and the DM is less cusped in the inner galaxy. 
Independently of the magnetic field used, the model is 
excluded if background uncertainties are of 1% (at 95% 
CL) or less, for DM masses below 40 GeV. Again, as 
commented above for the case of the effective approach, 
an equally valid option for the DM distribution, as can 
be the Einasto profile (not shown here) , can enhance the 
signal by a factor 3, allowing us to claim stronger con- 
clusions about the exclusion power of the synchrotron 
radiation. 

Finally, some words regarding the comparison with 
XENON100 bounds. It was shown in [78] that in general 
the scalar singlet model is disfavoured by XENON100 
data in the region of low (< Mw) masses, except 
around the region of the Higgs-pole, and also for m x < 
8GeV, where XENON100 starts loosing sensitivity. In 
this sense, synchrotron bounds are complementary to 
XENON100 bounds, as they are able to constrain the 
region of very low mass, as well as the region of large 
mass Mw), given some reasonable knowledge about 
the background. For illustration, we show in Fig. [TT] 
the exclusion of the coupling Cuss due to synchrotron 
data, assuming a value of 26/iG for the magnetic field 
in our ROI, and compare it directly with the exclu- 
sion coming from XENON100. We see that effectively 
for m x > Mw, synchrotron constraints are stronger 
than those of XENON100, assuming reasonable back- 
grounds uncertainties of 20% at 95%CL. Also, as ex- 
pected, synchrotron bounds (even in the conservative 
case of no-background) do better than XENON100, for 
masses m x < 12 GeV. 

One should notice that our analysis can also be applied 
to a vectorial dark matter through the Higgs portal [79] . 
Indeed, one can consider, to first approximation, a mas- 
sive vector as a scalar particle with three internal degrees 
of freedom. In this case the synchrotron flux should be a 
factor ~ 3 times greater than the one we obtained for the 
simplest singlet scalar extension. One thus deduces that 
the constraints we would obtain would be even stronger 
than the one we presented in the singlet scalar scenario. 



C. Extra U(l) 

Neutral gauge sectors with an additional dark U'(l) 
symmetry in addition to the Standard Model (SM) hy- 
percharge U(l)y and an associated Z' are among the 
best motivated extensions of the SM, and give the possi- 
bility that a dark matter candidate lies within this new 
gauge sector of the theory. Extra gauge symmetries are 
predicted in most Grand Unified Theories (GUTs) and 
appear systematically in string constructions. Larger 
groups than SU(5) or SO(10) allow the SM gauge group 
and U(l)' to be embedded into bigger GUT groups. 
Brane-world U'(l)s are special compared to GUT U'(l)'s 
because there is no reason for the SM particles to be 
charged under them; for a review of the phenomenology 
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FIG. 10. Top) Synchrotron flux at 45 MHz produced from the 
Higgs-portal model. For every mpjMi the coupling Ajjss is scanned 
while requiring that the relic abundance respects WMAP con- 
straint at 5cr. We presented the result for a magnetic field profile 
normalised to 26/iG (red, orange) and 10^tG (green, light green). 
Dashed line shows bound coming from data, assuming no back- 
ground; whereas dotted line shows the bound assuming a full back- 
ground with uncertainties of 5% (see text for details). Red and 
green bands are for NFW+MED+po = 0.43 GeV cm" 3 astro- 
physical setup, while light-green and orange correspond to ISO+la 
modcl+po = 0.3 GeV cm -3 . Bottom) Maximum background un- 
certainties allowed to exclude a point in the parameter space of the 
model, with the same conventions as above. 



of the extra U'(l)s generated in such scenarios see e.g. 
[80] . In such framework, the extra Z' gauge boson would 
act as a portal between the "dark world" (particles not 
charged under the SM gauge group) and the "visible" 
sector. 

Several papers considered that the key of the portal could 
be the gauge invariant kinetic mixing (5/2)Fy V F'^ |81j . 
One of the first models of dark matter from the hidden 
sector with a massive additional U'{\), mixing with the 
SM hypercharge through both mass and kinetic mixings, 
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FIG. 11. Constraints on the coupling Cf-^ss f° r the Higgs-portal 
model, as a function of DM mass m x . Different colours show the 
exclusion by synchrotron data, assuming: Red) Excluded by No- 
background choice, Blue) Excluded by uncertainties of 40%, Yel- 
low) Excluded by uncertainties of 20%, and Green) Excluded by 
uncertainties of less than 1%. Here we assume for the mag- 

netic field at GC, as well as NFW+MED set-up. Black solid line 
shows the exclusion by XENON100. 

can be found in [82] ■ The DM candidate tpo could be 
the lightest (and thus stable) particle of this secluded 
sector. Such a mixing has been justified in recent string 
constructions piHST] . but has also been studied within 
a model independent approach [83l [881 US] • For typical 
smoking gun signals in such models, like a monochro- 
matic gamma-ray line, see |90j . 

The matter content of any dark U(1)d extension of the 
SM can be decomposed into three families of particles: 

• The Visible sector is made of particles which 
are charged under the SM gauge group 5/7(3) x 
5/7(2) x /7y(l) but not charged under Up (1) (hence 
the "dark" denomination for this gauge group). 

• the dark sector is composed by the particles 
charged under /7d(1) but neutral with respect to 
the SM gauge symmetries. The dark matter (ipo) 
candidate is the lightest particle of the dark sector. 

• The Hybrid sector contains states with SM and 
Z7d(1) quantum numbers. These states are funda- 
mental because they act as a portal between the two 
previous sectors through the kinetic mixing they in- 
duce at loop order. 

From these considerations, it is easy to build the effective 
Lagrangian generated at one loop : 

l 1 8 

£ = -Csm — -^BiivB^ — -^X^X^ — -B^X^" 

i 3 



being the gauge field for the hypercharge, X^ the 
gauge field of Z7d(1) and ipi the particles from the hidden 
sector, VPj the particles from the hybrid sector, D„ = 
d,-i{q Y ~g Y B ll + q D ~g D X^gT a Wl), T a being the 5/7(2) 
generators, and 

with rrij and Mj being hybrid mass states [91] . 

Notice that the sum is on all the hybrid states, as they are 
the only ones which can contribute to the B^, X^ propa- 
gator. After diagonalising of the current eigenstates that 
makes the gauge kinetic terms of Eq. |29| ) diagonal and 
canonical, we can write after the 5/7(2)^ x U(l)y break- 
in£f] 

A,, = sin 9 w Wl + cos 9 W B^ (31) 
Z^ = cos0(cos0wW^ — smOwB^) — sincpX^ 
Z 1 = s\yl4>{cos9wW^ — sinflvK-Bju) + cos^A^ 

with, to first order in 8, 

a . 25sm0w 

cos <p = — - sin <p = — - 

V a 2 + 4<5 2 sin 2 9 W V a 2 + AS 2 sin 2 9 W 

a = 1 - Mf, /Ml - S 2 sin 2 9 W (32) 

± yj(l- M 2 ,/M 2 - S 2 sin 2 9 W ) 2 + AS 2 sin 2 9 W 

and + (-) sign if Mz' < (>)Mz- The kinetic mixing 
parameter 8 generates an effective coupling of SM states 
-0SM t° Z' , an d a coupling of ip to the SM Z boson 
which induces an interaction on nucleons. Developing 
the covariant derivative on SM and ipo fermions state, we 
computed the effective iPswPsmZ' and ipoipoZ couplings 
to first ordeiH in 8 and obtained 

£ = 9l>.9d(cos0 Z'^o^ipo + sin(p Z '">Po) ■ (33) 

We took qd9d = 3 trough our analysis, keeping in mind 
that our results stay completely general by a simple 
rescaling of the kinetic mixing 8. 

We show in Figure [T2] the synchrotron flux emitted at 
45 MHz for different values of Z' masses after a scan on 
(rriDM = rn^o',8) for points respecting the WMAP and 
electroweak precision tests (including Z width, p param- 
eter g — 2 and atomic parity violation constraint). Wc 
do this exercise for the astrophysical NFW+MED setup. 



16 Our notation for the gauge fields are (B^,X^) before the diag- 
onalization, (_B^,X^) after diagonalization and (Z M ,Z /Al ) after 
the electroweak breaking. 

17 One can find a detailed analysis of the spectrum and couplings 
of the model in the appendix of Ref. I92| . 
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We first observe a similar behaviour to that in the Higgs 
portal around two poles, Mz and Mz>- Indeed, the two 
channels giving a good relic density are the s— channel ex- 
change of the Z' bosons for tudm — Mz> /2. We have re- 
stricted the fluxes for each Mz> in windows to avoid jux- 
tapositions of the fluxes. However, when Mz> — 200 GeV 
the value of the flux begins to be quite weak for regions of 
other parameter space far from the pole region, due to the 
reduced factor in the amplitude square of the s— channel 
process ^ V>o ~> Z' -> SM SM \M\ 2 cx \jM%, which 
is almost a factor 1/6 compared to a 125 GeV Higgs ex- 
change, factor that we recover in the synchrotron fluxes. 
In this case one needs a precision better than 5% to be 
able to distinguish a dark matter signal if m,DM ^ 120 
GeV (see Figure 12 - bottom). 



One observes that the conclusions we have obtained for 
the Higgs portal are also quite valid for the Z' portal too, 
except that there are different fine-tuned pole regions for 
each mass of the Z' . The fluxes in both cases lie in a 
region of parameter space between 5% to 100% of the 
measured flux and could thus be probed in a near future. 



VI. CONCLUSIONS AND PROSPECTS 



In this work, we derived the constraints on synchrotron 
emission from the secondary products of dark matter an- 
nihilation in our Galaxy, based on the radio 45 MHz 
data. We expanded upon the current literature in sev- 
eral aspects: a) by using both a semi-analytical approach 
to model the particle propagation in the intergalactic 
medium [32J |33] , and a full numerical analysis based on 
the GALPR0P code [THJ[2D], benefiting from the strengths 
of both in particular aspects of the work, b) consider- 
ing intermediate Galactic latitudes, where the DM pro- 
files are more robustly constrained, and c) taking into 
account the large astrophysical uncertainties on the mag- 
netic fields. We have shown that synchrotron emissions 
can give very pertinent bounds on DM annihilation cross 
section, confirming previous results [T^l [T4THB] . 

The discussion of the astrophysical propagation setup, 
occupies a large fraction of this work. Together with a 
usual set of CR propagation parameters used previously 
to gauge uncertainties of the astrophysical conditions on 
the radio signals, (MIN, MED and MAX models of prop- 
agation [H]), we have considered several other sets of CR 
parameters, which are shown to be consistent with CR, 
gamma ray and/or synchrotron data [22ll24| . We have 
also explored the impact of the magnitude of the mag- 
netic field in our ROI and found that it is generally mild, 
ranging a factor of 3-4 for magnetic fields in the inner 
Galaxy in the range between 1-100 /iG. 

One of the purpose of our work is to put in perspective 
various indirect bounds on dark matter models, comple- 
mented with collider constraints. Although the scope of 
our results is broader, this was in part motivated by the 
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FIG. 12. Top) Scatter plot of prediction s fo r synchrotron flux 
coming from a Z'-portal model defined in l |29| . For every moMi 
the coupling 8 is scanned while requiring good relic abundance, for 
different Z' masses, m z i . A frequency of 45 MHz was used, assum- 
ing a magnetic field profile normalised to 26/iG (dark colours) and 
10/^G (light colours). Dashed line shows bound coming from data, 
assuming no background; whereas dotted line shows the bound 
assuming a full background with uncertainties of 5% at 95%CL. 
Bottom) Maximum background uncertainties allowed to be able to 
exclude a point in the parameter space of the model. 



fact that both colliders and indirect searches are sup- 
posed to put very relevant limits on relatively light dark 
matter candidates, ttidm 10 GeV, for which direct de- 
tection limits are weaker, or altogether inexistent. In this 
spirit, we have briefly reviewed (and applied to particle 
physics models in question, including forecasts) the con- 
straints from Fermi LAT observations of nearby dwarf 
spheroidal galaxies [T7] , as well as the limits set by CMB 
anisotropies [IS]. Both of these observations give strin- 
gent indirect limits on dark matter annihilation cross sec- 
tions and treat targets different than the one we focus on 
to derive the radio DM constraints, providing therefore 
independent probes of DM self-annihilation signatures. 
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We have applied our results to various DM particle 
models, starting with a generic approach based on effec- 
tive operators. There, we have shown that, for reason- 
able values of parameters, and for a conservatively chosen 
ROI, synchrotron searches for DM could be comparable 
to those of colliders, and sometimes even do better. In 
particular, in the case of effective couplings to leptons, we 
have shown that, even for very conservative setups, where 
we supposed that all the radio data are saturated by the 
synchrotron radiation produced by DM annihilation, the 
limits obtained are already competitive with those based 
on LEP measurements (photon+missing energy) . In this 
case too, radio constraints can even be better than the 
ones derived from CMB, depending on the uncertainty 
on the astrophysical modelling. Concretely, allowing that 
DM contribute to 5% the background signal could give 
the best limit obtained by current indirect detection sig- 
nals. For effective couplings to quarks, the synchrotron 
have a strong potential too, but due to huge uncertain- 
ties in the prediction of this signal, one cannot robustly 
claim that the most conservative synchrotron constraints 
can improve over the collider bounds. As expected, if 
DM couples to hadrons, dwarf constraint generally per- 
form better. 

The effective approach is powerful, but reaches its lim- 
its when resonance effects become relevant. To assess 
such effects, and for the own sake of a more microscopic 
approach to DM phenomenology, we have considered the 
constraints set by radio synchrotron radiation on two 
specific DM models: the so-called Higgs and Z' por- 
tals. This is motivated by the fact that these models 
are among the simplest extensions of the SM with DM 
candidates. Moreover, each of these models provide fully 
self-consistent UV completions of the effective operator 
approach, and thus provide complementary, albeit model 
dependent, information. We have shown that radio data 
may put severe constraints on these models, but at some 
price. Provided the uncertainty on the background could 
be assessed at the 5% level, and assuming a cuspy profile 
(NFW), and a large magnetic field, most of the Higgs 
portal parameter space is excluded, except near reso- 
nance (the Higgs pole), or close the threshold for W + W~ 
annihilation, which is impeded in the Galaxy (by con- 
struction both these effects are inexistent in an effective 
operator approach). We have also shown a very good 
complementarity between synchrotron bounds and the 
last XENON100 bounds when considering the exclusion 
of the parameter space of the model. For the case of the 
Z' portal, because of substantial annihilation into lep- 
tons, the constraints are stronger, or alternatively, the 
astrophysical setup may be more conservative (smaller 
magnetic field, less cuspy profile), but some control of 
the background is also required to exclude (most of) the 
parameter space. 

This work illustrates again, if necessary, that a multi- 
signal approach provides very complementary informa- 
tion on DM phenomenology. If on the long run, DM 



production at colliders is likely to give the strongest con- 
straints, one should keep in mind that missing energy 
may not be directly related to the actual DM that is 
supposedly present in our Galaxy. Our results present 
(part of) the state of art in confronting indirect searches, 
with a particular emphasis on synchrotron radio data, 
a very promising signal for dark matter, provided some 
control may be gained on the mundane, astrophysical 
background, a fascinating challenge for both future ob- 
servations and theoretical works. To pave the road, in 
the near future, PLANCK will be able to study Galactic 
emission in the frequency range where it is dominated 
by the dust emission, mapping with unprecedented pre- 
cision dust (and therefore indirectly gas) content of our 
Galaxy. Together with improvements in measurement of 
the charge cosmic ray spectra we will soon be getting 
from AMS-0i^1 and measurement of diffuse emission in 
gamma rays of such CR population, with the Fermi LAT, 
models of propagation and energy losses of CR are ex- 
pected to advance significantly over the next 5-10 year 
period. Finally, the future radio telescope facilities, like 
the LOFAR and SKA, will provide further leverage on 
the possible radio synchrotron signal from DM particles, 
in particular for lighter candidates. 
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Appendix A: Annihilation cross-sections (<jv), for 
different effective operators 



coming from the operators Os, Oy and Oa are given by: ^ I m 2 

°AV=-r^jJl--£ (A3) 



The expressions for the annihilation cross-sections 

48ttA 4 

" -n.ju-. (Al, * 24m / 




„ 8mt - 22m 2 m 2 f + 17mi „ 
24m 2 f + „ x \ L v 2 




cryv — i/l (A2) where m/ is the fermion mass and m x the DM mass, 

while v is the relative velocity of incoming particles, 
2 m 2 i g m 4 \ which at the moment of freeze-out it is assumed to be 

x f -^v 2 , v 2 = 0.24 0. 
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